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ABSTRACT 



Results from the X-ray spectral analysis of the ASCA PV phase observation of the 
Seyfert 1 galaxy IC 4329A are presented. We find that the 0.4 - 10 keV spectrum 
of IC 4329A is best described by the sum of a steep (F ~ 1.98) power-law spectrum 
passing through a warm absorber plus a strong reflection component and associated 
Fe K line, confirming recent results (Madcjski et al. 1995, Mushotsky et al. 1995). 
Further cold absorption in excess of the Galactic value and covering the entire source 
is also required by the data, consistent with the edge-on galactic disk and previous 
X-ray measurements. The effect of the warm absorber at soft X-ray energies is best 
parameterized by two absorption edges, one consistent with OVI, OVII or NVII, the 
other consistent with OVIII. A description of the soft excess in terms of blackbody 
emission, as observed in some other Seyfert 1 galaxies, is ruled out by the data. A large 
amount of reflection is detected in both the GIS and SIS detectors, at similar intensities. 
We find a strong correlation between the amount of reflection and the photon index, but 
argue that the best solution with the present data is that given by the best statistical 
fit. The model dependence of the Fe K line parameters is also discussed. Our best fit 
gives a slightly broad {a ~ 0.11 ± 0.08 keV) and redshifted (E ~ 6.20 ± 0.07 keV) Fe 
K hue, with equivalent width ~ 89 ± 33 eV. The presence of a weak Fe K hne with 
a strong reflection can be reconciled if one assumes iron underabundances or ionized 
reflection. We also have modeled the hne with a theoretical line profile produced by 
an accretion disk. This yields results in better agreement with the constraints obtained 
from the reflection component. 

Key-words: galaxies: individual (IC 4329A) - galaxies: Seyfert - X-rays: galaxies 
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1 Introduction 



X-ray observations of Seyfert 1 galaxies strongly support the hypothesis that optically 
thick (column density ^ 10^^ cm~^) matter surrounds the innermost regions of Active 
Galactic Nuclei (AGN) (Nandra & Pounds 1994). The detection of a bump in the con- 
tinuum emission above ~ 7 keV and an iron fluorescence line at ~ 6.4 keV in several 
Seyfert galaxies are compelling evidence that the intrinsic continuum is reprocessed by 
cold matter (Guilbert & Rees 1988, Lightman & White 1988). Recent ASCA observa- 
tions also confirm the existence of ionized matter (the so-called warm absorber) along 
the hue of sight of some Seyfert 1 galaxies (e.g., MCG-6-30-15: Fabian et al. 1994; 
NGC3227: Ptak et al. 1994; NGC3783: George, Turner & Netzer 1995). In these 
sources, the warm absorber alone can explain both the absorption features detected 
around 1 keV and the soft-excess emission. However, this seems not to be the rule for 
all Seyfert 1 galaxies (Matsuoka 1994). Some of them, hke NGC7469, do not show the 
presence of ionized absorption but require extra soft emission best characterized by a low 
temperature blackbody (Guainazzi et al. 1994). Others, like NGC4051 and Mkn766, 
require both the warm absorber and blackbody component (Mihara et al. 1994, Leighly 
et al. 1995, Guainazzi et al., in preparation). 

In light of these general considerations, the ASCA observation of the bright Seyfert 
1 galaxy IC 4329A (z=0.016, Wilson & Penston 1979) provides new insights on these 
issues. Prom the Ginga observation in the 1.7-30 keV energy band, IC 4329A is known 
to exhibit a strong reflection component which flattens the 1.7 - 30 keV photon index 
from an intrinsic value of 1.97 ±0.04 to the apparent value of 1.71 ±0.01 (Miyoshi et al. 
1988, Piro, Yamauchi & Matsuoka 1990, hereafter PYM90). The reflection component 
is accompanied by an iron fluoresence line at E ~ 6.2 keV, with an equivalent width 
(EW) of ~ 100 ± 20 eV. Evidence of day to day (~ 50%) variations in the amount of 
reflection relative to the direct continuum have also been reported (Fiorc ct al. 1992). In 
the ROSAT soft energy band (0.1-2.4 keV), the spectrum of IC 4329A is best described 
by an absorbed power law with a photon index F ~ 1.83 ± 0.4 plus an absorption edge 
around 0.7 keV having an optical depth of ~ 0.6 ± 0.2 consistent with absorption by 
OVI and OVII (Madejski et al. 1995). However, an alternative explanation of the soft 
excess in terms of a black-body or a steep power law is also allowed by the same data. 
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All fits derive a column density of a few 10^^ cm~^. This column density is largely in 
excess of the Galactic value A^Hgai — 4.55 x 10^° cm~^ (Elvis, Lockman & Wilkes 1989) 
and indicate absorption intrinsic to IC 4329A. This result is consistent with the fact 
that its host galaxy is seen edge-on (Petre et al. 1984). 

In the following, wc present the results from the X-ray observation of IC 4329A 
with ASCA (Tanaka, Inoue & Holt 1994). The data between 2 - 10 keV have been 
previously discussed in a letter by Mushotsky et al. (1995). We now report a more 
detailed analysis of the overall (0.4-10 keV) data. We show very strong evidence that 
the X-ray spectrum of IC 4329A is absorbed by highly ionized gas (warm absorber), 
characterized by the absorption edges of highly ionized oxygen ions (section 3.1). We 
also confirm the presence of a spectral hardening above ~ 6 - 7 keV, as well as the 
iron K-shell emission line, which are characteristics of a reflection spectrum (section 
3.2). Complex warm absorber models and accretion-disk emission line models are then 
discussed (section 4). We conclude that there is evidence of a warm absorber, a strong 
reflection component and a marginally broad iron fluoresence line in the X-ray spectrum 
of IC 4329A. 

2 The X-ray data 

The X-ray observation of IC 4329A reported here was carried out with ASCA on Au- 
gust 15th 1993 during the performance verification (PV) phase. Solid state imaging 
spectrometer (SIS) data were collected in 4-CCD 'Bright' and 'Faint' modes. Faint data 
were converted into Bright data. No echo corrections nor Dark Frame Error (DFE) cor- 
rections were applied (see Otani & Dotani 1994 for details on these effects). Standard 
selection criteria were used, the most relevant being elevation and bright Earth angles 
greater than 5° and 25°, respectively, and a magnetic cut-off rigidity greater than 8 
GeV/c. These were apphed to both CIS and SIS data. The resulting effective exposure 
times are ~ 29000 s for the SIS and ~ 32500 s for the CIS. About 55000 and 90000 
source plus background counts were collected from each SIS and CIS, respectively. The 
background count rates were about 1% of the SIS average source count rate (~ 1.7 
counts/s) and about 5% of the CIS average source count rate (~ 2.6 counts/s). 
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Only chip n.l of SISO and chip n.3 of SISl were used in the SIS analysis. This 
results in the loss of ~ 30 % of the total SIS source counts to the neighboring chips 
due to the broad wings of the mirror point response function. Source counts were 
extracted from circular regions of radius ~ 12 arcmin and ~ 6 arcmin for the gas 
imaging spectrometer (GIS) and SIS. respectively. Because the source is bright and fills 
the chip, SIS background was extracted from the SIS blank sky files from a same area as 
that used for the source. GIS background counts were accumulated from annular regions 
around the source, with an internal radius of ~ 14 arcmin and an external radius of ~ 
17 arcmin. Results with blank sky GIS backgrounds are consistent with the reported 
one to within ~ 5%. It has been reported from ROSAT observations (Madejski et al. 
1995) that two other X-ray sources (the Elliptical Galaxy IC4329 and an unidentified 
source) lie within ~ 13 arcmin of IC 4329A. In the following analysis, we will assume 
that their contribution is negligible since they become important only below ~ 0.5 keV, 
as pointed out by Madejski ct al. (1995). Neither source is spatially resolved with the 
SIS nor with the GIS, even when we considered only the soft energies (E ^ 0.5 keV) or 
only the hard energies (E ^ 5 keV). 

Although the source flux decreased by ~ 15% in a timescale of ~ 1 day in all 
detectors, a hardness ratio analysis reveals no significant spectral variation during the 
observation. Therefore, all photons have been accumulated together for the spectral 
analysis. 

3 Spectral analysis 

GIS2/3 and SISO/1 pulse height spectra were binned in order to have a number of 
counts/bin > 200 and ^ 100 in the energy ranges (0.7-10 keV) and (0.4-10 keV), respec- 
tively. Data preparation and spectral analysis have been performed using version l.Oh 
of the XSELECT package and version 8.5 of the XSPEC program (Arnaud et al. 1991). 
We made use of the response matrices gisv3_l.rmf (released in March 1994) for the 
GIS detectors and rspl.lalpha (released in June 1994) with Dark Frame Error (DFE) 
shifted by -2 ADU (- 7 eV) and -5 ADU (- 18 eV) for the SISO and SISl detectors, 
respectively. 
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For the GIS data, we first modelled the overall (0.7 - 10 keV) spectrum with a single 
absorbed power law plus iron emission line. The column density, A^h, was free to vary. 
This model gives an acceptable (x^g^/d.o.f. ~ 1.33/564) description of the spectrum 
with r ~ 1.79 ± 0.02, a column density Nn (3.8 ± 0.2) x lO^^cm-^ and a very broad 
{a ~ 0.53±0.24 keV) iron emission line at E ~ 6.19 ±0.15 keV with EW ~ 233 ± 80 eV 
(see Table lA). With these values, the unabsorbed flux in the 2-10 keV energy band 
is ~ 1.01 X 10~^° erg cm~^ s~^, which corresponds to a luminosity L2_iokeV ~ 1-07 x 
10^^ erg s~^. Although the continuum shape is approximately consistent with previous 
measurements with the same model (PYM90, Madejski et al. 1995), the iron hne appears 
much broader and a factor of ~ 2 stronger than previous values (PYM90, Nandra & 
Pounds 1994). However, this difference can be easily understood if we consider that, as 
shown later on, the SIS data strongly requires the addition of a reflection component 
to the continuum emission (section 3.2). Because of its better resolution, the SIS can 
indeed better separate the Fe K line from the reficction component. When the reflection 
is included in the model, the Fc K line parameters obtained from the GIS data are 
considerably reduced (Table IB) and become consistent with Ginga and SIS results. A 
detailed description of the reflection parameters given in Table IB, as well as the full 
justification for the use of this model, is presented in section 3.2, in the hght of the 
results obtained from the SIS spectrum. 

For the SIS data, a single absorbed power law model with Ny{ ~ 4 x 10^^ cm~^ and 
r ~ 1.75 gives an unacceptable fit, with x^g^/d.o.f. ~ 2.21/551. The residuals to this 
fit between 0.4 - 10 keV are shown in Figure 1 and are indicative of the complexity of 
the X-ray spectrum of IC 4329A. At soft energies, the residuals show the presence of 
excess emission at E < 0.8 keV and absorption features around 1 keV. These indicate 
either that the underlying soft continuum emission is more complex than a single power 
law or that the absorption is more complex than a uniform, cold absorption. At higher 
energies, we find an Fe K emission hne at E ~ 6.2 keV but also systematic positive 
residuals starting from E > 5 keV. The presence of this "hard tail" will appear more 
evident (section 3.1.2.) after fitting adequately the soft energy part of the spectrum. 

Because of the complexity of the SIS spectrum, we split the following analysis in 
two parts: (i) fit of the spectrum between 0.4 - 5 keV in order to focus only into the 
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characteristics of the soft energy band, (ii) fit of the overall (0.4-10 keV) spectrum. 

3.1 Soft energy band 
3.1.1 Thermal component 

One possibility to explain the measured soft excess would be to assume a two compo- 
nent continuum emission instead of a single power law. A steep, ultra soft excess may 
be associated with emission from the inner parts of an accretion disk (Wilkes & Elvis 
1987, Turner and Pounds 1989). As a first approximation, this soft emission can be 
parameterized by a low temperature (T ^ 10^ K, kT ^ 100 eV) blackbody. A model 
consisting of a blackbody plus absorbed power law was fitted to the 0.4 - 5 keV SIS 
spectrum (SISO and SlSl are fitted simultaneously). The resulting blackbody temper- 
ature of ~ 9.3 X 10^ K (kT ~ 80 eV), photon index of ~ 1.96 and column density of 
~ 5.8 X 10^^ cm~^ are consistent with previous ROSAT results (Madejski et al. 1995), 
but the residuals shown in Figure 2a present strong and systematic deviations around 
E ~ 0.6 keV and absorption features around 0.9 keV. These clearly indicate that the 
blackbody model cannot, alone, explain the measured soft excess. The further addition 
of an absorption edge at E ~ 0.76 keV with r ~ 0.46 improves the statistics of the fit to 
X^g^/d.o.f. ~ 1.22/481. In that case, the normalization of the blackbody is reduced by 
a factor of ~ 5 and some systematic structures are still present in the residuals around 
0.9 keV. These imply that also this possibility is unlikely. 

Alternatively, we modelled the soft excess with a Raymond-Smith plasma model 
(Raymond & Smith 1986) absorbed by the Galactic column only. This model has been 
proposed in the past as a possible explanation for the soft excess observed in some Seyfert 
1 galaxies (Turner et al. 1991). The best result is obtained allowing the abundances to 
be free. This gives a plasma temperature of ~ 2.1 x 10^ K (kT ~ 183 eV) and improves 
the fit to x^g^/d.o.f. ~ 1.28/479. However, the best-fit abundances are ~ 1/20 the 
solar value. Such a required low abundance is obvious from the absence of any blend 
of emission lines around 0.8-0.9 keV. Moreover, some systematic features are present in 
the residuals (Figure 2b). Thus we also exclude this alternative as a possible description 
of the soft X-ray spectrum of IC 4329A. 
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3.1.2 Absorption edges 

Another possibility to explain the soft excess is an ionized or "warm" absorber. The 
ionized absorber model assumes that light elements (O, Ne, Mg, Si, Fe) arc photoionized 
and therefore predicts extra emission due to the reduction of the opacity at soft energies. 
This model also predicts extra absorption around 1 keV, mainly due to highly ionized 
oxygen (Netzer 1993). As clear from the residuals shown in Figure 1, both these features 
are detected in the spectrum of IC 4329A. These strongly support, a priori, the ionized 
absorber model. 

As a first approximation of a warm absorber, we added a single absorption edge to 
the absorbed power law model. Given the good statistics of the present data (~ 50000 
counts/SIS) and the ASCA spectral resolution (FWHM ~ 80 eV at 1 keV), a single 
absorption edge is a rather crude approximation of a warm absorber, but it allows a 
direct comparison with previous measurements. The inclusion of an absorption edge 
in the absorbed power law model improves the fit to x^g^/d.o.f. ~ 1.30/479 and gives 
Eedge — 0.76 ± 0.01 keV and Tedge — 0.84 ± 0.04. These values are consistent with the 
ROSAT results (Madejski et al. 1995). However, inspection of the residuals in Figure 
2c clearly suggests further complexity. Indeed, the addition of a second absorption 
edge improves the statistics to xled/^-'^-^- — 1-15/477, which is significant at more 
than 99.9% level (Ax^ ~ 76). As evident from the residuals in Figure 2d, the two 
edges plus absorbed power law model provides the best description of the soft X-ray 
spectrum. Best-fit parameters are reported in Table 2 and confidence contours for the 
edge energies versus optical depths are shown in Figure 3. Assuming a conservative 
value of ~ 10 eV {la) for the SIS systematic uncertainties, the edge at E ~ 0.70 ±0.02 
keV is consistent with absorption by OVI(671 eV), OVII(739 eV) and NVII(667 eV), 
where the energies quoted in parenthesis are the edge energies in the rest frame (Lotz 
1968). The best agreement is with OVII. The second edge at E ~ 0.84 ± 0.02 keV 
is consistent with 0V111(871 eV), only. Inspection of the residuals in Figure 5 shows 
the presence of a further edge-like feature at ~ 1.25 keV, consistent with NeIX(1196 
eV), NeX(1362 eV) and Fe L edges of Fe ionization stages around FeXVIII. Inclusion 
of a third absorption edge into the model is, however, significant only at a ~ 68 % 
confidence level. Given its low statistical significance, we do not consider it any further 
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in the following. It is interesting to note that both characteristic energies obtained for 
the OVII and OVIII edges are redshifted by approximately the same amount (~ 25 
eV, in the observer frame). We checked that this redshift is not related to the Charge 
Transfer Inefficiency (CTI) of the SIS (see Otani & Dotani 1994 for details). 

Because IC 4329 A is a fairly bright point source, wc may expect any systematic 
problems due to remaining calibration uncertainties to be more visible in IC 4329A 
than in other, weaker sources. We point out that some features are present between 
~ 1.8 - 2.5 keV. These are likely to be due to calibrations uncertainties in this energy 
range (Si K-edge and Au M-edge of the detector). In this work, we will assume that 
these features are not of cosmic origin, and are therefore of no physical interest. 

Also, even after applying our best- fit edge model (above), there are still residuals 
in the shape of an emission line around 0.6 keV and an absorption feature around 0.5 
keV. In particular, the emission line can be fitted with a narrow {a = 0) gaussian 
line at 0.58 ± 0.05 keV with EW ~ 30 ± 25 eV. The inclusion of this emission 
line corresponds to a Ax^ ~ 4, which is significant at more than 95%. Its energy is 
consistent with emission from OVII and/or OVIII. The existence of similar emission 
lines is predicted in warm absorber models (Netzer 1993) and has been confirmed in the 
Seyfert 1 galaxy NGC3783 (George, Turner & Netzer 1995). However, it is also known 
from other ASCA observations that features around 0.5-0.6 kcV arc very sensitive to 
calibration uncertainties (due to the O K-edge at ~ 0.54 keV), background substraction 
and selection criteria (Miura et al. 1995). Therefore, we have carefully checked to 
estabhsh the believability of this soft emission fine. 

When different background regions are used and different selection criteria (hke 
BR_EARTH > 50 for ex.) are applied, both the absorption and soft emission line 
features remain unchanged. Using the "gain" command in XSPEC, we (linearly) shifted 
the energies of SISO and SISl response matrices by ± 2 ADU (~ 7 eV), i.e. within typical 
systematic errors of ~ 10 eV (la) for the DFE "zero-level" (Otani & Dotani 1994). After 
fitting the data again with the shifted matrices, the absorption feature around 0.5 keV 
appears smoothed out. It is therefore probably due to calibration uncertainties. While 
the soft emission line changes in intensity it remains statistically significant in every 
case. However, the emission fine is significant only in the SISl data. Also, it appears to 
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go away when only the night data are considered although, in this case, the statistics 
are considerably reduced. We conclude that with the present data, the existence of this 
emission line is doubtful. 

3.2 Overall spectrum 

We now include also the data from E > 5 keV and proceed to the analysis of the 
spectrum in the broad (0.4 - 10 keV) energy band. One striking result is that we detect 
a clear hardening of the spectrum even down to 5-6 keV, likely due to the presence 
of a reflection component, as previously found by Mushotsky et al. (1995). Figure 4 
illustrates the existence of the hard tail and of a Fe K emission line around 6.2 keV. To 
produce this Figure, the spectrum was fitted between 0.4 - 5 keV with the two edges 
plus power law model (section 3.1.2., Table 2) and the data above 5 keV were added 
back in. The data were rebinned to more than 500 counts/bin in order to better show 
the presence of the hard tail. The systematic positive values of the residuals, before 
and after the FeK emission line, are a clear indication that the continuum is increasing 
at higher energies. This is also confirmed by the huge (cr ~ 1.17 keV, EW ~ 760 eV) 
values obtained for the Fe K line (E ~ 6.25 keV) if we include only a single gaussian 
line in the model, which tries to account for the line plus the hard tail. We underline 
that in this well as in all the following spectral fits, all the parameters of the 

two edges plus power law in the model are free. Given the importance of the residuals 
above 7 keV, we carried a number of tests to confirm the reahty of the hard tail. The 
use of a background choosed from the outer region of the relevant chip in the SIS had no 
significant effect on either the overall spectrum or the strenght of the hard tail. The use 
of different extraction radii for the source (from r=4 arcmin to r=9 arcmin) and more 
conservative selection criteria (like BR_EARTH>50 to avoid light leakage) gave similar 
results as well. 

In a recent ASCA long observation of the Seyfert 1 galaxy MCG-6-30-15, a very 
broad iron line has been observed (Tanaka et al. 1995). It has been well modelled by 
the sum of a broad (E ~ 5.6 keV, a ~ 0.54 keV, EW ~ 150 eV) and narrow (E ~ 6.5 
keV, a ~ 0.20 keV, EW ~ 150 eV) gaussian lines. These have been interpreted as the 
signature of the "double- horned" line expected from an accretion disk (e.g., Fabian et 
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al. 1989). If a double-gaussian parameterization is used for IC 4329A instead of a single 
line, one gaussian converges to a very broad {a ~ 1.2 keV, EW ~ 713 eV) line centered 
at ~ 6.2 keV and the other one to a narrow {a — 0, EW ~ 53 eV) at the same energy. 
Since the data require an unusually broad gaussian that is non-physical and inconsistent 
with any previous X-ray measurement, we prefer to interpret this result as meaning that 
the high energy part of the spectrum requires an additional continuum component. 

The presence of a Fe K emission line and a hard tail is usually interpreted in terms 
of a model in which the primary continuum emission is reprocessed by optically thick 
matter, i.e. a reflection component (Lightman & White 1988, Nandra & Pounds 1994). 
Therefore, we added a reflection component in the two edges plus power law model. 
A gaussian emission line is also included in the fit. The refiection used is the "plrefl" 
model in XSPEC and is calculated following Lightman & White (1988). The reflecting 
matter is assumed to be thick, cold and with cosmic abundances (Morrison and Mc 
Gammon 1983). The inclination angle i has been flxed to 0, i.e. the reflecting plane is 
seen face-on, and the covering factor VL/2ti has been flxed to 1. The source is assumed 
to emit isotropically and the spectrum is integrated over all viewing angles. The relative 
normalization R ( = Arefl/Api) between the reflected component (A^e/;) and the direct 
component (Ap^) is a free parameter. This implies that the normalization of the reflection 
component is allowed to be independent of the normalization of the power law. The 
parameter R, alone, includes the uncertainties concerning the geometry of the reflector 
(i.e., the covering fraction), the possibility that the source is anisotropic (e.g., Ghisellini 
et al. 1991) and time-lag effects. This imphes, for example, that a best-flt value of R 
~ 1 can be explained if the reprocessing material covers ~ 27r of the source. If R > 
2, then the result cannot be explained only in terms of a covering factor increase, and 
either the source must be anisotropic or there are time delay effects. 

The photon index has been left free to vary. The best-fltting parameters obtained 
with a reflection flt are given in Table 3. The unfolded spectrum and residuals are shown 
in Figure 5. The measured photon index of 1.98 ± 0.05, A^h ^ (3.48 ± 0.12) x 10^^ 
cm-2 and Fe K line parameters (E ~ 6.20 ± 0.07 keV, a ~ 0.11 ± 0.08 keV and EW 
~ 89 ± 33 eV) are consistent with previous Ginga results (PYM90, Fiore et al. 1992). 
Confldence contours in the Fe K parameter space cr - E are well constrained (Figure 6). 
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These show that the Fe K Une is consistent with redshifted cold iron at more than 90% 
confidence level, and is marginally broad at an ~ 80% confidence level. As discussed 
later on, these characteristics of the iron line suggest that we are dealing with a disk 
geometry (Fabian et al. 1989, Mushotsky et al. 1995). 

As already mentioned in section 3, the GIS spectral analysis is completely consistent 
with that derived above for the SIS data. In fact, the inclusion of a reflection component, 
assuming an inclination i — and covering factor Jl/27r = 1 as for the SIS spectrum, 
improves the GIS fit at a confidence level > 99% (Ax^ ~ 21). Gomparison of the 
GIS best-fit parameters (Table IB) with those obtained from the SIS data (Table 3), 
shows very good agreement from both detectors and allow us to be fairly confident in 
the proposed best-fitting model. The major discrepancies are in the A^h and photon 
index best-fit values which are respectively higher and fiatter in the GIS than in the 
SIS spectrum. However, these differences are likely attributed to the fact that the GIS 
is not sensitive enough at soft energies to detect the absorption edges. As they are 
not fitted explicitely in the GIS spectrum, they likely enhance the absorption value and 
therefore contribute to fiatten the GIS spectrum. This hypothesis is also supported by 
the fact that the fit of the spectrum with all 4 detectors (SIS -|- GIS) gives results almost 
identical to that found with the SIS only. 

3.3 On the model-dependence of photon index, reflection nor- 
mahzation and Fe K hne parameters 

With the above best statistical fit of IC 4329A (Table 3), we found that there is a strong 
(linear) correlation between the photon index and the amount of refiection (Figure 7). 
This relation should be borne in mind when interpreting the results of this analysis. A 
similar effect has been discussed for the analysis of MGG-2-58-22 (Weaver et al. 1995). 
While a high best-fit value of R may indeed indicate that the true intrinsic photon index 
is steep, it also depends strongly on fitting the correct amount of refiection, as well as 
the shape of the refiected spectrum modelled. There is no guarantee that the current 
refiection model is entirely correct and we may biasing the fits with a slightly incorrect 
shape or normalization of the reflection. Different analysis procedures may, therefore, 
lead to slightly different results. For example, flxing the photon index to its 90% upper 



12 



limit (r ~ 1.92) obtained from the fit in the 0.4 - 5 keV energy band (Table 2), the 
amount of reflection converges to R ~ 2.1. However, flxing the photon index from the 
flt at soft energies is not necessarily correct in this case because we may also bias our fit 
with an improper parameterization of the warm absorber. Alternatively, fixing R = 1, 
a value appropriate for a face-on slab subtending 27r solid angle at the X-ray source, 
the fit gives F ~ 1.89. However, also in this case, there is no reason to fix R. The best 
solution of the present data is therefore that obtained with F and R free to vary (Table 
3). 

Because of the correlation between photon index and the amount of reflection in the 
model we investigated by how much this could affect the parameters obtained for the Fe 
K line. For this, we calculated the 90% confidence contours of a versus energy for Fe K 
with R varying from zero (no reflection) to 3.35 (the best-flt value). These contours are 
shown in Figure 8 for R=0, 0.3, 1, 1.5, 2.1 and 3.35. The corresponding best-fit values 
of the Fe K line width are 1.17, 0.89, 0.21, 0.15, 0.13 and 0.11 keV, respectively. The ^ 
values are 641, 640, 639, 634, 630 and 626 for 547 d.o.f. We see that the energy of the 
Fe K line depends only slightly on the exact amount of refiection. These contours are 
indicative of 3 different "solutions" for the line width: (i) assuming httle or no reflection 
(R ^ 0.9) results in a very broad Fe K line with characteristic width of ^ 0.8 keV and 
EW ^ 380 eV, (ii) assuming an intermediate contribution from reflection (R ~ 1 - 2) 
results in two local minima, one near o ~ 0.6 - 0.8 keV and the other near cr ~ 0.1 - 
0.2 keV, and (iii) assuming a strong reflection component (R from ~ 2 to the best-flt 
value ~ 3.3) results in one minimum at a ~ 0.11 keV with well constrained contours. 
Fitting simultaneously the 2 SIS and 2 GIS detectors yielded very similar results. These 
contours, as well as the best-fit values, are similar to those obtained by Mushotsky et al. 
(1995) although their fit yields a slighlty larger Fe K width than inferred from our best- 
fit (case (iii)). The difference may be a consequence of the slightly different continuum 
modelling, probably caused by the different energy band used by the authors. 

To summarize, we find that although the photon index and reflection normalization 
are strongly coupled together, the most likely description of the spectrum of IC 4329A is 
given by its best statistical flt (Table 3). We have also shown that despite this intrinsic 
uncertainty, the Fe K parameters depend only little on the exact amount of reflection. 
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provided R is higher than ~ 2, as required by the best- fit model. With these caveats 
borne in mind, we can discuss and try to interpret these results physically. 

4 Discussion 

Above results show that the best interpretation of the soft X-ray spectrum of IC 4329A 
is in terms of photoionized absorber. Therefore, we have also fitted the SIS data between 
0.4 and 5 keV with a warm absorber model, also including a uniform cold absorber. The 
model has been produced using the photoionization code CLOUDY (Ferland 1991) and 
using the assumptions as in Fabian et al. (1994). Free parameters are the ionization 
parameter ^ = h/nB? (erg s~^ cm), the "warm" column density Ny^ (cm~^), and the 
cold column density A^h (cm~^). The photon index has been fixed to 2.0, a value close 
to the SIS best- fit value. The fit yields an acceptable fit with x^g^/d.o.f. = 1.23/480. 
Although it is statistically better than with the blackbody model, it is worse than with 
the two absorption edges. We attribute this difference to the fact that the model used 
is based on simple basic assumptions which may not reflect the reality, namely: a thin 
geometry for the photoionized gas, solar abundances, constant temperature, constant 
density and a source luminosity of 10^^ erg s~^. Further, the model does not include the 
emission lines expected from a more reahstic warm absorber (Netzer 1993). We obtain 
log ^ ~ 1.12 ± 0.04, log A'w ^ 21.70 ± 0.02 and log A'h ^ 21.49 ± 0.01. The fitted 
neutral column density is again about one order of magnitude higher than the Galactic 
value. 

Interesting is to note that with the 2 edges parameterization of the warm absorber, 
both edges best-fit energies are redshifted by approximately 25 eV. If real, this redshift 
correspond to a Doppler velocity of about 10000 km s~^. Although this result depends 
strongly on the correct modelling of the warm absorber (e.g., the redshift of the lower 
edge may be interpreted as a blend of OVI and OVII edges) and on the precise cahbration 
of the energy scale, it may be the signature of infalling matter or largescale turbulance 
associated to the warm absorber region. Firm conclusions on this point require more 
investigations, but we point out that edge energies shifted compare to their expected 
energy have been reported also for the observations of MCG-6-30-15 and NGC4051 
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(Otani 1995). 

As suggested by Krolik & Kriss (1995), resonance line scattering may significantly 
contribute to the total opacity of the warm gas. It is worth noting therefore that it 
is a potential confusing factor when interpreting the energy and optical depth of the 
absorption edges. For example, the presence of blended resonance oxygen lines may 
cause a misidentification of OVIII as OVII due to the apparent energy shift in the onset 
of the OVIII edge. Alternatively, it could contribute to shift the apparent energies of the 
absorption edges. Unfortunately, quantitative considerations on resonance absorption 
being strongly model-dependent (e.g. on the covering factor and velocity of the warm 
gas), and given the limited spectral resolution of the SIS, we cannot exclude the presence 
of resonance absorption lines in the present data. 

A striking alternative to the warm (plus cold) absorber is absorption by cold matter 
with abundances free to vary. IC 4329A is known to exhibit absorption of a few 10^^ 
cm~^. This absorption is consistent with the fact that its host galaxy is seen edge-on, 
and so we examine whether the absorption edges seen in the ASCA data can be due to 
neutral interstellar matter. Since the K-edge energies of OVII, OVIII and NeX are close 
to the L-edge energy of Fel, K-edge energies of Nel and Mgl, respectively, it is possible 
that abundances larger than the cosmic value can mimic the absorption edges normally 
interpreted as the presence of a warm absorber. In fact, a fit of the soft energy band 
(0.4 - 5 keV) with a cold absorber with variable abundances gives an equally good fit as 
that obtained with our previous parameterization using two warm edges (section 3.1.2). 
However, extremely high abundances (~ 5 - 15 times the cosmic value) of all metals 
are required. Also, known cahbration problems around 0.5 - 0.6 keV (see section 3.1.2) 
make the results from this fit doubtful. A deeper investigation using these and other 
data is currently being conducted and will be presented in a future paper (Hayashida 
et al., in preparation). 

The X-ray spectrum of IC 4329A also requires a remarkably strong reflection com- 
ponent. In a recent re-analysis of HEAO — 1 (2-40 keV) data, a similar amount (R 
~ 2.3 -ir 4.2) of refiection component has been detected in IC 4329A (Weaver, Arnaud, 
& Mushotzky 1995). Similarly, we find in the ASCA spectra ~ 2 - 3 times more re- 
flection relative to the continuum than in the Ginga observation (PYM90, Fiore et al. 
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1992, Nandra & Pounds 1994). The present data also suggest a slight, but significant, 
increase in the absolute amount of reflection. Although we might recall that the ASCA 
best-fit spectra, as well as the HE AO — 1 and Ginga spectra, exhibit a strong coupling 
between photon index and amount of reflection (see section 3.3), these results indicate 
the presence of a strong and variable reflection component. It should be noted that a 
time-resolved analysis of the Ginga observation also suggested day to day variations in 
the amount of reflection component in IC 4329A (Fiore et al. 1992). 

If the geometry of AGN consists on a source illuminating hard X-rays above and 
below a soft X-ray emitting disk (e.g. Haardt, Maraschi & Ghisellini 1994), and these 
hard X-rays are produced by Inverse Compton scattering of the soft X-rays, then the 
radiation field within the hard X-ray emitting region is likely to be strongly anisotropic 
(Ghisellini et al. 1991). The reflecting matter may therefore see more radiation than 
that emitted towards the observer, by up to a factor of ~ 5, as a direct consequence 
of the electrons losing more energy on head-on collisions with the soft X-ray photons 
than on collisions from behind. In that strong reflection component becomes 

physically plausible. However, anisotropy alone is unlikely to explain the increased con- 
tribution of the reflected component compared to that of the direct continuum because 
it requires a change in the anisotropic factor. In other words, anisotropy can explain a 
high value of R but not a variation of R. Instead, a time lag effect could, alone, explain 
all of the observational constraints. Though IC 4329A has never exhibited strong, rapid 
variabihty, it has varied up to a factor of ~ 2.8 (from a 2-10 keV flux of ~ 5 x 10~^^ erg 
cm~^ s~^ to ~ 14 X 10~^^ erg cm~^ s~^) in a timescale of ~ 5-6 years, with a doubling 
timescale down to ~ 1 year, in its ~ 20 years history of X-ray observations. A time lag 
could therefore cause strong amount of reflection. Our results can indeed be understood 
in this context if the continuum flux decreased by some amount prior to the ASCA 
observation. Specifically, the flux would have to have decreased by a factor of ~ 2 from 
a time when it was higher than that measured by Ginga, and the reflection component 
could not have had time to respond to this decrease. This can simultaneously explain a 
high value of R, a variation of R, and a higher absolute amount of reflection as compared 
to Ginga. 

However, another difficulty with the above results is that strong reflection should be 
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accompanied, in most cases, by a strong iron line. We note that the reflected spectrum 
used in our fitting procedure (section 3.2) is averaged over all viewing angles and has 
accordingly a shghtly different shape than a face-on reflected spectrum. The best-fit 
amount of reflection R ~ 3.3 corresponds therefore to a (~ 1.33 times) lower value R ~ 
2.5 for a facc-on view (Ghisellini, Haardt & Matt 1994, Magdziarz & Zdziarski 1995). 
As a rough estimate, assuming an Fe K line EW of ^ 140 eV in the case R=l for a 
cold face-on disk and an incident power law spectrum with F ~ 2.0, we expect an EW 
of ~ 300 eV. The best-fit amount of reflection requires therefore ~ 3 times more line 
than observed (EW ~ 89 ± 33 eV). This difference may be reconciled by reducing the 
abundance of iron, or by increasing the abundance of the light metals (lighter than iron) 
(Reynolds, Fabian & Inoue 1995). Alternatively, if the Fe K line arises in an ionized 
disk, its EW can be very low even in the presence of strong reflection because of either 
resonant trapping of line photons or complete ionization of iron (Ross & Fabian 1993, 
Zycki & Czerny 1994). Compare to the neutral case, the ionized disk has lower opacity 
at lower energies. This might help to explain why the reflection shows up at such low 
(~ 5-6 keV) energies in this source. However, we show below that the above constraints 
can be considerably reduced because the EW obtained with an accretion disk line model 
is significantly higher. 

The parameters of the Fe K line obtained from the SIS spectral analysis arc consistent 
with the line being redshifted by ~ 100 eV and slightly broad {a ~ 110 eV). These 
support the idea that a fit of the hne with a disk line model (Fabian et al. 1989) would 
be more appropriate than a gaussian fit. A disk line model takes into account the effects 
of doppler-broadening and gravitational and transverse redshifts in order to calculate 
the characteristic line profile. The resulting line profile is highly asymmetric, and often 
"double-horned" . A remarkable example of the applicability of this model to Seyfert 1 
galaxies has come from ASCA long observations of MCG-6-30-15 (Tanaka et al. 1995). 

We attempted a fit of the SIS data with the "diskline" model available in the XSPEC 
program. This model calculates the intensity and profile of the line from a flat, optically 
thick accretion disk as calculated in Fabian et al. (1989). The details of the line depend 
on the inclination i of the disk, the inner Rj^ and outer Rout disk radii, and on the 
hne emissivity a (assumed to vary as R~"). Following Mushotsky et al. (1995), we fix 
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Rj„ = 3 Rg and Kout — 1000 Rg, where Rg is the Schwarzschild radius of the black hole 
(Rg ~ 3 X 10^M/Mq cm, with M the mass of the black hole). The inclination of the disk, 
i, is also fixed at (face-on) to be consistent with our assumption of a face-on reflection 
component. The line energy is fixed to 6.3 which correspond to the redshifted energy 
of a fluorescence Fe K line. Only the normalization of the line and its emissivity a are 
free to vary. The fit yields a ~ 2.30 ± 0.20 and an EW of 188 ± 60 eV. If i is allowed 
to vary, we obtain an upper limit of i ^ 25 degrees which confirms a posteriori that the 
assumption of a disk face-on is substantially correct. In this case, the EW is enhanced 
to 213tll eV. These results are characteristic of a disk observed face-on, in which most 
of the fluorescence Fe K line is emitted within a few tens of Schwarzschild radii. As 
stated above, with the use of a disk line model, the Fe K EW is in rough agreement 
with the value expected from the reflection component. 

Strictly speaking, if a disk line model is included, we should also take into account 
relativistic effects on the (continuum) reflection component in order to have a self- 
consistent model. This model is not available at the moment. This problem will be 
explored in a later paper. However, relativistic corrections should produce merely an 
energy shift and/or broadening of the reflection component. Therefore, we expect these 
shouldn't significantly affect the above results. 

It should be noted however that the velocity width (FWHM ~ 12000 ±8500 km s"^) 
of the Fe K line is still consistent with an average derived FWHM of ~ 4070 ± 30 km 
s~^ for the strong Balmer lines (H^, H^ and H-y) which dominate the optical spectrum of 
IC 4329A. This suggests that the Fe K hue might, alternatively, be produced from the 
broad line region. One-sided beaming of the continuum emission and/or obscuration 
may be necessary, however, to explain the apparent redshift of the line. 

5 Conclusions 

To summarise, the present observation has revealed the following. 

(i) The apparent soft X-ray spectrum of IC 4329A is complex with a predominant 
soft-excess (compared to a simple absorbed power law) at E < 0.8 keV and deep absorp- 
tion structures between ~ 0.8 - 1.2 keV. The unprecedent high (> 10 in every channel) 
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signal to noise ratio of this observation allows us to discard, directly from the residuals, 
a blackbody emission model as an alternative explanation for the soft excess. The soft 
X-ray spectrum is instead best described by two absorption edges at E ~ 0.70 keV and E 
~ 0.84 keV, consistent with absorption by OVI, OVII or NVII, and OVIII, respectively. 
This complex absorption is clearly the signature of a warm absorber, confirming and 
extending previous results obtained from ROSAT observations (Madejski et al. 1995). 

(ii) Confirmation of the existence of a "hard tail" at higher energies, characteristic 
of a refiection component (Mushotsky et al. 1995). The hard tail is clearly detected 
in both GIS and SIS, at similar intensities. When the refiection is included, we find 
that the data prefer a stronger absolute amount of reflection and stronger contribution 
relative to the direct continuum than previously observed. A time lag effect can account 
for all observational constraints. 

(iii) Because of the presence of the refiection component, the Fe K line width 
and equivalent width are model-dependent. These are, however, only weakly model- 
dependent provided R is higher than ~ 2.1 (section 3.3). With the proposed best-fit 
model (with R ~ 3.3 ± 1.1), we find evidence that the Fe K fine is redshifted (by ~ 100 
eV) and slightly broad {a ~ 0.11 ±0.08 keV), as expected if the emission fine is emitted 
by an accretion disk. The weak (EW ~ 89 ± 33 eV) Fe K (gaussian) fine requires iron 
undcrabundance or/and ionized matter to be consistent with strong reflection. A flt 
with a diskline model which includes the wings of the line can, however, account for 
most of the discrepancy as well. 
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TABLE 1 - GIS2+3 Data (0.7 - 10 keV) 
A. Results of Power-Law + Fe K Line Fits 
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B. Results of Power Law + Fe K Line + Reflection Fits 
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Note: Intervals are at 90% confidence for one interesting parameter. 
" Normalization in 10~^ photons cm""^ s"^ keV"^ at 1 keV. 
R = Arefi./Api (see section 3.2 for details) 



TABLE 2 - SISO+1 Data (0.4 - 5 keV) 
Results of Two Edges + Power Law Fits 
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Note: Intervals are at 90% confidence for one interesting parameter. 
" Normalization in 10~^ photons cm"^ s"^ keV"^ at 1 keV. 
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TABLE 3 - SISO+1 Data (0.4 - 10 keV) 
Results of Two Edges + Power Law + Fe K Line + Reflection Fits 
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Note: Intervals are at 90% confidence for one interesting parameter. 
" Normalization in 10~^ photons cm~^ s~^ keV~^ at 1 keV. 
R = Arefl./Api (see section 3.2 for details) 
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Figure captions 

Figure 1: Residuals of the SIS spectrum between 0.4 - 10 keV, for a model consisting 
of a single absorbed power law with F ~ 1.75. This figure illustrates the soft excess, 
absorption features around 1 keV, Fe K hne at ~ 6.2 keV and hard tail starting at E > 
5 keV (see also Fig. 4). 

Figure 2: Residuals of the SIS spectrum between 0.4 - 5 keV for a model consisting 
of: (a) absorbed blackbody plus power law, (6) absorbed power law plus Raymond 
Smith spectrum covered by only the Galactic column density, (c) absorption edge plus 
absorbed power law and {d) two absorption edges plus absorbed power law. 

Figure 3: Contour plots illustrating the (observed) energy and optical depth 68, 90 

and 99 per cent confidence levels of the two absorption edges at (a) E ~ 0.7 keV and 
(6) E ~ 0.84 keV. The redshifted energies of the absorption edges have been marked in 
the figures. 

Figure 4: Residuals of the SIS spectrum between 0.4 - 10 keV, for a model consisting 
of a two absorption edges plus an absorbed power law fitted between 0.4 - 5 keV. The 
data have been rebinned with more than 500 counts/bin to better show the systematic 
increase of the residuals above 5 keV. 

Figure 5: SISO+1 unfolded spectrum, best-fitting power law plus two absorption edges 
plus Fe K line plus reflection model (upper panel). Residuals in the form of the ratio 
(data/model) are shown in the lower panel. 

Figure 6: Fe K width and energy 68, 90 and 99 per cent confidence contour levels for 
the best fit model, as given in Table 3. 

Figure 7: Contours at 68, 90 and 99 per cent confidence of amount of reflection (R) vs. 
photon index (F). They illustrate the strong correlation between the two parameters. 

Figure 8: Contours at 90 per cent confldence of amount of hne width, cr, vs. observed 
line energy calculated for various amounts of reflection. The best-flt values are marked 
with crosses. 
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TABLE 3 - SISO+1 Data (0.4 - 10 keV) 
Results ol Two Edges + Power Law + Ee K Line + Reflection Eits 
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" Normalization in 10 ^ photons cm ^ s ^ keV ^ at 1 keV. 
^ R = Ai.efi./Api (see section 3.2 for details) 



